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ABSTRACT 

Calcium-rich supernovae (Ca-rich SNe) are peculiar low-luminosity SNe Ib with rela¬ 
tively strong Ca spectral lines at ^2 months after peak brightness. This class also has 
an extended projected offset distribution, with several members of the class offset from 
their host galaxies by 30 - 150 kpc. There is no indication of any stellar population 
at the SN positions. Using a sample of 13 Ca-rich SNe, we present kinematic evidence 
that the progenitors of Ca-rich SNe originate near the centers of their host galaxies 
and are kicked to the locations of the SN explosions. Specifically, SNe with small pro¬ 
jected offsets have large line-of-sight velocity shifts as determined by nebular lines, 
while those with large projected offsets have no significant velocity shifts. Therefore, 
the velocity shifts must not be primarily the result of the SN explosion. There is an 
excess of SNe with blueshifted velocity shifts within two isophotal radii (5/6 SNe), 
indicating that the SNe are moving away from their host galaxies and redshifted SNe 
on the far sides of their galaxies are selectively missed in SN surveys. Additionally, 
nearly every Ca-rich SN is hosted by a galaxy with indications of a recent merger 
and/or is in a dense environment. At least 6-7 host galaxies also host an AGN, a 
relatively high fraction, again linking the nuclear region of these galaxies to explosions 
occurring tens of kpc away. We propose a progenitor model which fits all current data: 
The progenitor system for a Ca-rich SN is a double white dwarf (WD) system where 
at least one WD has a significant He abundance. This system, through an interac¬ 
tion with a super-massive black hole (SMBH) is ejected from its host galaxy and the 
binary is hardened, significantly reducing the merger time. After 10 - 100 Myr (on 
average), the system explodes with a large physical offset. The rate for such events is 
significantly enhanced for galaxies which have undergone recent mergers, potentially 
making Ca-rich SNe new probes of both the galaxy merger rate and (binary) SMBH 
population. 

Key words: supernovae—general, supernovae—individual (SN 2000ds, SN 2001co, 
SN 2003H, SN 2003dg, SN 2003dr, SN 2005E, SN 2005cz, SN 2007ke, 
SN 2010et, SN 2012hn, PTFOOdav, PTFllbij, PTFllkmb), galaxies—individual 
(2MASS J22465295-b2138221, CGCG 170-011, IC 3956, NGC 1032, NGG 1129, 
NGG 2207, NGG 2272, NGG 2768, NGC 4589, NGC 5559, NGC 5714, NGC 7265, 
UGC 6934) 


1 INTRODUCTION 

For the past century, systematic supernova (SN) searches 
have discovered thousands of SNe. Nearly all of these SNe 
fall into 3 classes: Type la. Type II, and Type Ib/c. How¬ 
ever, over the last decade with the implementation of large 
SN searches, we have begun to discover many astrophys- 
ical transients that do not fall into the well-delineated 
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classes mentioned above. It turns out that there were 
about a dozen new classes of “exotic” or “peculiar” tran¬ 
sients lurking in the s hadows. These cl asses include lu¬ 
minous SNe I ln (e.g.. Smith et al.l 120071 ) . T ype I super- 
luminous SNe lO uimby et al.|[201ll ) , kilonova e (IBerg er et al.l 


I2OI3I : iTanvir et al. 20131). SNe lax dFolev et alT 20131) . and 
SN 2006bt-like SNe ( Folev et al.ll201ol ). 

One class, the SN 2005E-like SNe llPerets et al.l[201ol . 
hereafter. iPlOl ). also known as “Calcium-rich SNe” (or sim- 
ply “Ca-rich SNe”; iFilippenko et all l2003l : iKasliwal et al.l 
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|2012| . hereafter, lK12h . are particularly interesting. These 
SNe are technically of Type Ib, having distinct He lines 
in their spectra near maximum brightness. However, they 
are l e ss lum inous and faster fading than normal SNe Ib 
iPil : IKI 2 I) and are of ten found in early-type galaxies 
(PIC; iLvman et al .11201^ . Additionally, several members of 
the class, including SN 2005E have large projected offsets 
( '>10 kpc and up to at least 150 kpc) from their h ost galaxies 
IpIOI : IK12| : lYuan et al.ll201 iLvrtian et al.l[20l3 : this work). 

Relative to other SNe Ib, Ca-rich SNe also have a dis¬ 
tinct spectroscopic evolution. After about 2 months, the SNe 
show strong forbidden lines, indicating that the ejecta are 
at that point mostly optically thin. This is a much faster 
transition than for most SNe Ib. This “nebular” spectrum 
is also distinct, showing extremely strong [Call] AA7291, 
7324 emission relative to that of [Q I] A A6300, 6363, giving 
the name to the class (iFilippenko et alll2003ll . Modeling of 
the ejecta also indicates that the m ass-fr action of Ca in the 
ejec ta is also extremely high 1^1/3: [PlOll . 

IPIOI examined the rate of Ca-r ich SNe in the Lick Ob - 
servatory Supernova Search (LOSS: [Pilippenko et al.|l200lll. 
Hnding that they occur at 7 ± 5% the rate of SNe la. ml 
found a rate >2.3% that of SNe la for the Palomar Transient 
Factory (PTF) sample. Therefore, Ca-rich SNe are some¬ 
what common events and cannot come from extremely rare 
progenitor scenarios. _ 

Since most C a-rich SNe have early-type hosts ilS; 
iLvman et al. I [mi), a massive star origin is unlikely. More¬ 
over, there is no indication of star formation at the posi¬ 
tion of the SNe in deep pre- or post-explosio n imaging (IPlOl : 
IPerets et al.ll20fll : lK12l : iLvman et ^l2014ll . 

Because of the large offset distribution, it has been 
suggested that C a-rich SNe occur in d warf galaxies (|K12; 
[Yuan et "ni 120131 ) or globular clusters (lYuan et al.l 1 20131) : 
however, the deep limits on any stellar light at the positio n 
of the SNe rule out such possibilities (ILvman et al.l 1^141) . 
These observations require that the progenitors of Ca-rich 
SNe be born elsewhe re and travel a sign ificant distance to 
where the y explode dLvman et al.ll2014l) . With this deter¬ 
mination, |Lxma^^eF3!n2^^~suggested that Ca-rich SN 
progenitors were neutron star (NS) - white dwarf (WD) bi¬ 
nary systems which are kicked by the SN that created the 
NS and then undergo a merger after traveling far from their 
birth site. 

In this manuscript, we present a kinematic and host- 
galaxy study of Ca-rich SNe. By examining the velocity shift 
distribution, the projected offset distribution, and the angle 
offset distribution, as well as the host-galaxy properties of 
the sample, we constrain the progenitor systems and origin 
of Ca-rich SNe. In Sections [5] and [31 we dehne our SN sample 
and discuss their host-galaxy properties. In Sections [4] and 
(5) we measure the projected offset and velocity shift distri¬ 
butions for our sample. In Sectionwe analyze these data. 
We present a basic progenitor model in Section[3 We discuss 
our results and summarize our conclusions in Section [8] 


2 SUPERNOVA SAMPLE 

Our sample consists of all known (and pub lished) Ca-rich 
SNe. We begin with the ILvman et al] (l2014l) sample of Ca- 
rich SNe. This sample includes 12 SNe originally reported 


in other works dPlOl: Kawabata et al. I I2OIC1I : ISullivan et al.l 
I2OIII : IKI2I : IValenti et al. 2014h . We include PTF09dav in 
this analysis, although at maximum light, it had much 
lower ejecta velocities than other members (ISullivan et al.l 

B and uniquely displays hydrogen emission at late-times 


To this sample, we add PTFllkmb. PTFllkmb was 
discovered on 24.24 August 2011 (all times are UT) by PTF 
l|Gal-Yam et al.ll201ll) . A Keck/LRIS spectrum obtained 28 
August 2011 wa s used to classify PTFllkmb as a SN Ib 
dCal-Yam et al.ll201ll ). We use this spectrum as well as a 
Keck/LRIS spectrum from 26.37 November 2011 to classify 
PTFllkmb as a Ca-rich SN fAppendixIXll. 

Basic host-galaxy information was obtained through 
the NASA/IPAC Extragalactic Database (NED) and is pre¬ 
sented in Table \M\ Four SNe have somewhat ambiguous 
hosts, and we discuss each below. 

SN 2003H occurred between two merging galaxies, 
NGC 2207 and IC 2163. Their recession velocities differ by 
24±25 km s“^, and this difference does not affect our results. 
SN 2003H is offset by 8.73 and 5.75 kpc from NGG 2207 and 
IG 2163, respectively, and the difference in offset does not 
affect any results. 

SN 2007ke occurred in a cluster environment (AWM 7) 
near the brightest member, NGC 1129. It is also near a group 
member, MCG-l-07-07-003. The difference in their recession 
velocities is 227 ± 23 km s“^, with NGC 1129 having the 
larger redshift. This difference does not affect any results, 
but is noteworthy. The SN was offset by 16.71 and 8.20 kpc 
from NGC 1129 and MCG-l-07-07-003, respectively. Again, 
this difference does not affect any results. Given the relative 
sizes of these galaxies and large offset from MGG-l-07-07- 
003, NGC 1129 is a reasonable choice for the host galaxy. 
Nonetheless, given the cluster environment, it is possible 
that SN 2007ke originated from a faint cluster member or 
was formed in the intracluster medium. 


SN 2010et occurred 37.6 - 69.2 kpc from three galax¬ 
ies at redshifts consistent with the SN redshift. The closest 
galaxy, CGCG 170-011 is also relatively large. Using the 
NED reported major and minor 2MASS Kg isophotal axes 
(with a reference value of Kg = 20.0 mag arcsec“^) and po¬ 
sition angle, SN 2010et was offset 5.5 isophotal radii from 
CGCG 170-011. SN 2010et is 69.2 kpc and 16.5 isophotal 
radii from the other nearby, large galaxy, GGCG 170-010. A 
third galaxy, SDSS J171650.20-1-313234.4, is 46.9 kpc and is 
too small and faint for a 2MASS isophotal measurement. 
Based on offset, the most likely host is CGCG 170-011, 
and we assume this for the rest of the analysis. 112 de¬ 
termined that no source exists at the position of SN 2010et 
with Mh < —12.1 mag, but this does not completely rule 
out a possible dwarf galaxy as the host of SN 2012et. 

Finally, PTFllkmb is far from any galaxy. PTFllkmb 
has a redshift of roughly g = 0.017 as determined by SN fea¬ 
tures (iGal-Yam et al.ll201ll . and conhrmed in our analysis). 
The closest galaxy listed in NED is 4.2 arcmin away and has 
no redshift information. However, if it is at the distance of 
PTFllkmb, it would be offset by ~80 kpc. The three clos¬ 
est galaxies with redshifts are 2MASX J22224094-I-3613514 
(offset by 4.5 arcmin, 83.0 kpc, and 36.6 isophotal radii; 
cz = 4500 ± 25 km s“^), UGG 12007 (offset by 6.5 arcmin, 
127.3 kpc, and 15.0 isophotal radii; cz = 4829±33 km s“^), 
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Figure 1. DSS images of Ca-rich SN host galaxies, with each centered in individual panels. Each panel is 8^ X 8^ except for PTFllkmb, 
which is 16' X 16'. All panels have the same scale and are aligned with North up and East to the left. Each SN position is marked with 
a blue circle. An ellipse at twice the isophotal radius is shown in red. 


and NGC 7265 (offset by 7.3 arcmin, 150.0 kpc, and 10.2 
isophotal radii; cz = 5083 ± 26 km s“^). 

While NGC 7265 is very far physically from PTFllkmb, 
it is the closest in terms of isophotal radii and has the most 
similar redshift to that derived from the SN. We therefore 
adopt NGC 7265 as the host of PTFllkmb for this analysis. 

All spectra ( except for that of PTFllkmb) were previ¬ 
ously publ i shed (iKaw abata et al.ll201(il : IPIOI : ISullivan et al.l 
l20li] : IK12I : IValenti et al.l l2 0l4l. We obtained some of these 
data through WISERep llYaron fc Gal-Yatnll2012ll . 

Almost all SNe in our sample have only one late-time 
spectrum. For the few SNe with multiple late-time spectra, 
we primarily analyze the latest high-quality spectrum avail¬ 
able. The spectra typically had phases of ~2 months after 


maximum brightness (e.g., Hi). At these phases, and for 
all spectra examined here, Ca-rich SNe have strong [Call] 
AA7291, 7324 emission, especially relative to that of [Ol] 
AA6300, 6363. 


3 HOST-GALAXY PROPERTIES 

3.1 Evidence for Recent Host-galaxy Mergers 

In Figure [T] we present Digitized Sky Survey (DSS) images 
of the host galaxies of the Ca-rich SN sample. The host 
galaxi es of Ca-rich SNe are c learly atypical as previously 
noted (IPld : iLvman et al.|[201^ ; the sample is highly skewed 
to early-type galaxies. 
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However, previous studies had not noted the strong ev¬ 
idence for recent host-galaxy mergers. In fact, nearly every 
Ca-rich SN host galaxy shows some indication of a recent 
merger and/or is in a very dense environment where the 
likelihood of recent galactic mergers is much larger than in 
the field. 

Of the 13 host galaxies, one is clearly inter¬ 
acting (NGC 2207) and one i s a disturbed sp i- 
ral (2MASS J22465295-f2138221; ISullivan et ^ l201lD . 


NGC 1129 h as a disturbed morp hology indicative of a re¬ 
cent merger dPeletier et al.lE99Cll l. There are also four SO 
galaxies (IC 3956, NGC 1032, NGC 2272, and NGC 7265), 
whic h are often though t to be the result of recent mergers 


Moore et al 


Moore et al 


199S ) or which may have been “harassed” 


19981 ). From these tracers alone, 7/13 Ca- 


(e.g., 

(e.g., 

rich host galaxies have some indication of a recent merger. 

Additionally, at least 8/13 host galaxies are in a group 
or cluster environment. This is an exceedingly high fraction, 
even for elliptical galaxies. Moreover, 7/13 host galaxies are 
either the brightest group galaxy (BGG) or brightest cluster 
galaxy (BCG). As these galaxies tend to sit at or near the 
center of the group/cluster potential, they are more likely 
to have had recent mergers than typical cluster members. 

The relative rate of Ca-rich SNe to SNe la is roughly 
7% ([pT 3). Using the relative fraction of “recent merger” host 
galaxies for the two groups, the Ca-rich SNe would have a 
relative rate of only 2 ±2% that of SNe la in galaxies with no 
indications of recent mergers and 11 ± 8% that of the SN la 
rate in galaxies with some indication of a recent merger. 

In total, 11/13 Ca-rich SN host galaxies either have 
some evidence for a ongoing/recent merger or are in dense 
envir onments. Using the LOSS SN la sample (iLeaman et al.l 
[2Qi3), we find that 52% of SN la host galaxies have simi¬ 
lar host properties. Using binomial statistics, there is only a 
1.6% chance that SNe la and Ca-rich SNe come from similar 
host populations in terms of their “recent merger” proper¬ 
ties. Since Ca-rich SN and SN la host galaxies have indistin- 
guishable morphology an d star-formation rate distributions 
l|P10l : iLvman et ^120131 '). the merger/environment is likely 
a significant factor in the relative rates. The galaxies in rich 
environments or that have had recent mergers have an en¬ 
hanced rate of Ca-rich SNe relative to that of SNe la. 


3.2 A High Occurrence of AGN 


In addition to the different indications of merger activ¬ 
ity, several host galaxies also host an active galactic nu- 
cleus (AGN). Pr eviously identified AGN include NGC 2207 
jKjmfman et al.l l2012h. NGC 276 8 JV eron-Cettv fc V eronI 
l2006h . NGC 4589 Jho et al. 1 19971 : iNaear et al.l 20051 ). and 
IC 3956 fe.g.. iLiu et al .1 120111) . Additionally, NGC 1032 has 
a radio-to-far infr ared flux ratio consistent with an AGN 
dPrake et al.ll200^ . We also examined the Sloan Digital Sky 
Survey (SDSS) emission li ne quantities for the host galax¬ 
ies observed by SDSS fe.g.. iTremonti et al.ll200^ . For these 
five galaxies, we identify GGCG 170-011 as a likely AGN, 
confirm that IC 3956 has a “composite” spectrum indicative 
of AGN activity, and identify NGG 5559 as having a likely 
composite spectrum (Figure [2]). 

At least 6-7 (depending on the classification of 
NGC 1032) of the 13 Ca-rich host galaxies also host an 
AGN. This is a large fraction considering that this is an 



log ([N II]/Ha) 

Figure 2. Bald win-Phillips-Terlevich (BPT) diagram 
llBaldwin et al.l Il98ll) for the Ca-rich SN host galaxies with 
SDSS spectra. The curves delineating the star-forming region 
(below the dashed line), the AGN region (above the solid line) 
and the “comp osite” region (betwe en th e dashed and solid lines ) 
are taken from lKewleY _^ alj (l200ll) an d lUauffmann et all (l2003l ) 
as implemented in iKewlev et al.l ll200fil ). 


incomplete survey of Ca-rich host galaxies. Although the to¬ 
tal AGN fraction may be as high as -^43% when the lowes t 
luminosity galaxies and AGN are included le.g.. Ih3|2008I ). 
this fraction can be lower for specific galaxy samples. For in¬ 
stance, an SDSS spectroscopic sample of >120,000 galaxies 
found that 18% of all galaxies hosted an AGN as determined 
from line diagnostics llKauffmann et aLll2003l) . For the host 
galaxies with SDSS spectra, 3/5 galaxies also host an AGN. 
Using binomial statistics for this subsample with the SDSS 
occurrence rate, there is a 4.7% chance that the Ca-rich host 
galaxies are typical. If we expand this to the larger, incom¬ 
plete sample, we find a <0.4-2.0% chance (depending on if 6 
or 7 galaxies host an AGN). However, we note that if we use 
a 43% occurrence rate, these differences are not significant 
for any subsample. For such a rate, ^9/13 galaxies must 
host an AGN for a significant difference. 

There is likely a connection between the nuclear regions 
of these galaxies and the SNe offset up to 150 kpc in projec¬ 
tion. Since the typical duty cycle of an AGN is ~10^ years, 
the information from the nucleus to SN location would have 
to travel at a very high speed of >1000 kms”*^. 


4 GALACTIC OFFSETS 

Ga-rich SNe have been found very far from their host galax¬ 
ies fe.g.. |P10l : lK12l) . However, this is not exclusively the case. 
For instance, SN 2003dg was offset 4.1" from its host galaxy 
corresponding to a projected distance of only 1.7 kpc. 

Using NED values for the angular size scale and offsets 
between the SNe and host-galaxy, we determine a projected 
offset for each SN. The typical uncertainty in the offsets is 
^.1", corresponding to projected offset uncertainties of only 
0.01 - 0.07 kpc. All SN properties, including those relative 
to their host galaxies, are presented in Table 
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In Figure [T] we present Digitized Sky Survey (DSS) im¬ 
ages of the host galaxies of the Ca-rich SN sample. In each 
image, we display an ellipse corresponding to two isophotal 
radii for each galaxy and the SN position. 

To determine the isophotal radial offset for each SN, we 
use the 2MASS Ks isophotal elliptical parameters (with a 
reference value of Ks = 20.0 mag arcsec“^) for each galaxy. 
For each SN position, we determine a multiplicative factor 
such that an ellipse with the same axis ratio, position angle, 
and center will intersect with the SN position. This mnlti- 
plicative factor is the number of isophotal radii at which the 
SN is offset. 

There are many sources for isophotal parameters, but 
2MASS provides a large and homogeneous source. However, 
these parameters may have somewhat large errors for par¬ 
ticular galaxies. The most suspect value is the position angle 
for NGC 2207. The 2MASS position angle is 70°, while other 
sources have values of 97° and 141°. It appears that these 
differences are caused by the different position angles of the 
bulge/bar and the spiral arms. Nonetheless, NGC 2207 has 
an axis ratio of 0.68, so the exact angle has only a marginal 
affect on the measured isophotal radial offset for SN 2003H. 

Some galaxies are highly inclined, which is important 
for understanding the SN offsets. For instance, SN 2003dr 
is offset only 2.7 kpc (in projection) from its host galaxy 
nucleus, but it is offset nearly perpendicular to its nearly 
edge-on host galaxy (minor axis relative to major axis of 
0.16). Thus SN 2003dr is offset by 2.9 isophotal radii from 
the nucleus of its host galaxy. 

4.1 Galaxy Targeted SN Search Efficiency 

Some Ca-rich SNe have extremely large angular offsets. For 
instance, SN 2005E and PTFllkmb are offset 2.3' and 7.3' 
from their host galaxies, respectively. Galaxy-targeted SN 
searches, such as LOSS, could potentially miss several Ga- 
rich SNe because they are offset beyond the field-of-view 
(FOV) of the camera. This is particularly worrisome since 
PTF, which runs an untargeted search and has a large FOV 
camera, has only discovered (and announced) Ca-rich SNe 
that have a projected distance of at least 34 kpc and up to 
150 kpc. 

To better understand the efhciency of the LOSS search 
at finding Ca-rich SNe, we perform a simple calculation. 
The KAIT camera has a 7.8' x 7.8' FOV. By assuming that 
a targeted galaxy is centered in the middle of the FOV, we 
can determine the fraction of projected offsets observed for 
a galaxy at a given distance. 

Figure [3] displays the results from this calculation. For 
very nearby galaxies, where the scale is 0.1 kpc/", KAIT 
observes 10, 50, and 100% of the area within 83, 37, and 
23 kpc, respectively. Since many Ca-rich SN e have peak ab¬ 
solute magnitudes of M ~ —16 JPlOl : [k 1^ . KAIT will be 
magnitude limited for Ca-rich SNe beyond -^80 Mpc. At 
this distance, the scale is ~0.4 kpc/". For this scale, KAIT 
observes 10, 50, and 100% of the area within 333, 149, and 
93 kpc, respectively. A typical Ca-rich host galaxy within 
the KAIT sample may have >^0.3 kpc/", where KAIT ob¬ 
serves 10, 50, and 100% of the area within 250, 112, and 
70 kpc, respectively. 

The farthest offset in the Ca-rich SN sample is 
PTFllkmb with an offset of 150 kpc. Of the remaining PTF- 



Projected Offset (kpc) 


Figure 3. Fraction of the region within a certain projected off¬ 
set from a particular galaxy observed by KAIT. The dotted, 
long-dashed, short-dashed, and solid lines correspond to scales 
of 0.1, 0.2, 0.3, and 0.4 kpc/", respectively. The arrows represent, 
from smaller to larger offsets, the projected offsets of SN 2005E, 
PTF llbij, and PTFllkmb, respectively. 


discovered Ca-rich SNe, the median offset is 37.6 kpc, and 
therefore KAIT should have been able to discover such SNe 
in nearly all of their targ eted galaxies. SN 2 005E, which 
was discovered by LOSS llGraham et al.l [20051 1. was offset 
by 24.3 kpc. In fact, all Ca-rich SNe, with the exception 
of PTFllkmb, are offset less than 3.9', meaning that they 
would land within the KAIT FOV. 

Although KAIT will not be able to discover SNe with 
extremely large offsets, especially for the closest galaxies, it 
should be able to sample the bulk of the population discov¬ 
ered by PTF. Although there are likely systematic effects 
related to humans identifying a transient far offset from a 
host galaxy as a true SN, there is no technical reason for a 
lack of Ca-rich SNe with large offsets in the LOSS sample. 

Since LOSS is not clearly missing Ca-rich SNe with 
large offsets, it is reasonable to assume that PTF is missing 
some Ca-rich SNe with small offsets (or not reporting these 
objects as Ca-rich SNe). This may be caused by the diffi¬ 
culty of detecting faint sources on top of small galaxies, the 
difhculty of obtaining high-quality spectra of faint sources 
on top of relatively bright galaxies, or a preference to ob¬ 
tain multi-epoch, high-quality spectra of SNe Ib with large 
offsets. 

PTF not detecting or not classifying some Ca-rich SNe 
projected on top of galaxies may also explain the “low” 
rate for Ca -rich SNe relative to the KAIT rate measurement 

i|pTol :li^ . 


5 VELOCITY SHIETS 

Forbidden lines, which are produced by optically thin ma¬ 
terial, are extremely useful for understanding the velocity 
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distribution of SN ejecta along our line of sight. In partic¬ 
ular, forbidden emission lines provide the opportunity to 
measure the bulk velocity offset of the SN ejecta along our 
line of sight. 

For Ca-rich SNe, there are two optical forbidden line 
complexes which may be studied, [O I] AA6300, 6363 and 
[Call] AA7291, 7324. As a defining property of this class, 
the SNe have much stronger [Call] lines than [Olj. As a 
result, the majority of our sample has weak or undetectable 
[O I] lines. However, the [Ca II] lines are sufficiently strong 
for proper kinematic measurements in all spectra used for 
this analysis. 


5.1 Velocity Measurement Method 

For a single emission line, there are three typical measure¬ 
ments used to determine the bulk velocity offset: the peak 
of emission, the center of a Gaussian fit to the emission 
profile, and the emission-weighted velocity. If a line profile 
is Gaussian, the same velocity would be measured from all 
three methods. However, if there is significant skewness to 
the profile or multiple peaks, the methods could differ. 

For [Call] AA7291, 7324 in SN spectra, the doublet is 
not typically resolved, and this is the case for Ca-rich SNe 
as well. However, the lines are separated enough such that 
a single Gaussian is generally not a good description of the 
data. Moreover, the line profiles are generally moderately 
skewed, making even the sum of two Gaussian profiles (off¬ 
set to match the wavelength offset of the doublet) a poor 
description of the data. 

As an example. Figure |4] shows the [Call] AA7291, 7324 
doublet for SN 2005E. For all line profiles, we subtract the 
underlying continuum by linearly interpolating across the 
line profile. The shape of the profile is complicated and far 
from Gaussian. The profile peaks blueward of the nominal 
central wavelength with a large drop in emission around zero 
velocity. Although this ha s been an indicati on of dust for¬ 
mation in SN ejecta (e.g., ISmith et al.ll200^ . the long tail 
of redward emission (extending perhaps 3000 km s”*^ fur¬ 
ther to the red than the blue) suggests that the line profile 
cannot be explained simply by dust reddening. Instead, the 
profile is likely indicative of a complex velocity structure for 
the SN ejecta. 

Clearly a single Gaussian is not an appropriate descrip¬ 
tion of the [Gall] profile for SN 2005E. We have also at¬ 
tempted to fit the profile with two Gaussians where their 
rest-frame wavelength offsets correspond to the difference in 
wavelengths for the doublets, and their widths and heights 
are constrained to be the same. These constraints should re¬ 
sult in a perfect match to the data if the emitting material 
were described by a simple velocity distribution. However, 
SN 2005E has a more complex velocity distribution. Because 
of the poor Gaussian fits for many spectra in our sample, 
we do not use this method for our analysis. 

A second possibility is to use the peak of emission. Since 
most spectra contain noise at the level that would affect re¬ 
sults (or possibly small scale structure), we smooth the spec¬ 
tra before determining the peak. For SN 2005E, the velocity 
of the peak of emission is offset from the best-fit Gaussian 
velocity by 1100 km s~^. The peak of the emission is a par¬ 
ticularly poor tracer of the entire SN ejecta since it lacks any 
information about the shape of the velocity distribution, can 
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Figure 4. Continuum-subtracted [Call] AA7291, 7324 profile for 
SN 2005E (black curve). The spectrum is also reflected across 
zero velocity (dotted grey curve) to demonstrate the skewness in 
the tails of the feature. A best-fit double Gaussian profile with 
rest wavelengths corresponding to those of the [Call] doublet is 
shown as the red dashed curve. The velocity shifts corresponding 
to the peak of the smoothed line profile, the emission-weighted 
center of the profile, and the shift of the double-Gaussian profile 
are marked with small vertical blue, gold, and red lines. 


be influenced by small density perturbations in the ejecta, 
and is strongly affected by noise in the spectra. We therefore 
do not use this method of our analysis. 

The final possibility is the emission-weighted velocity. 
Relative to this velocity, half of the emission is blueshifted 
(and half redshifted). This method uses all information from 
the line profile, is not affected by the two emission lines, 
is not significantly affected by noise, and is not affected 
by skewness. Moreover, this velocity has a specific physical 
meaning. If the ejecta are completely optically thin and each 
Ga atom in the ejecta has an equal probability of emitting 
a photon at these wavelengths, this velocity will correspond 
to the systemic line-of-sight velocity. As this is clearly the 
best velocity measurement for determining any line-of-sight 
velocity shifts, we use this method for our analysis. 


5.2 Systematic Uncertainties 

While the method chosen to measure of the velocity shift 
is relatively robust, we investigate possible systematic un- 
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certainties, which dominate the uncertainty of the velocity 
measurement. 

To determine the uncertainty, we performed a Monte 
Carlo simulation for each spectrum. For each spectrum, we 
smoothed the [Ca II] line profile and measured the emission- 
weighted velocity. Then using noise properties matched to 
the data, we simulated 1000 spectra for each spectrum. We 
measured the emission-weighted velocity of each simulated 
spectrum, randomly changing the regions used to determine 
the continuum. There was no bias in these measurements, 
with most shifts from the noise-free spectrum being one or 
two pixels. This demonstrates the robustness of this method 
to noisy spectra and with different choices for the contin¬ 
uum. We use the standard deviation of the shifts from the 
noise-free measurement as the uncertainty in the emission- 
weighted velocity. 

There is an additional potential systematic uncertainty 
related to the continuum subtraction. Our method assumes 
a linear underlying continuum. However, if the continuum is 
dominated by a galaxy spectrum, a linear continuum is not 
necessarily appropriate. An example of likely host-galaxy 
contamination for a SN spectrum is shown in Figure[Sl where 
we show the spectrum of SN 2003dg. 

We examine the effect of subtracting a galaxy contin¬ 
uum rather than a linear continuum. For this task, we choose 
an elliptical template spectrum, which appears to be an ex¬ 
cellent match to the continuum for SN 2003dg, and a red¬ 
dened Sc template spectrum since its host galaxy is classified 
as Scd. Using the elliptical and Sc template spectra result 
in velocity shifts that are systematically 240 and 80 kms“^ 
bluer than the linear continuum resulting in a velocity shift 
of —1360 and —1200 kms“^, respectively, compared to our 
nominal value of —1120 kms“^. 

We have examined all spectra, and only two other 
spectra potentially have host-galaxy contamination. For 
SN 2001co, the continuum flux is low enough that the choice 
of continuum does not result in any difference in velocity 
shift. For SN 2003dr, both the elliptical and Sc templates 
result in a systematic blueshift of 80 kms“^. 

While it is possible that our choice of continuum is bias¬ 
ing the measured velocity shifts for SNe 2003dg and 2003dr, 
we choose to use the results from having a linear contin¬ 
uum. This choice is conservative given our results (the other 
choices strengthen our findings). It also provides an easier 
path to reproducing the met hodol ogy with future datasets. 

Following the method of IPld (and references therein), 
we also attempted to subtract a “photospheric” SN contin¬ 
uum from each spectrum assuming that the underlying con¬ 
tinuum is caused by photospheric SN emission. We matched 
our spectra to those of other SNe Ib with similar spectral fea¬ 
tures away from the nebular lines. We then used these spec¬ 
tra as the underlying continuum. Doing so did not change 
the measured velocity in any case. For i nstan ce, performing 
the same continuum subtraction used by|pT3 for SN 2005E, 
the measured velocity offset shifted by <1 pixel. 

This is not surprising given the relative strength of the 
nebular lines to the continuum. In the case of SN 2005E, 
which has one of the strongest continua, the continuum is 
only ~15% the flux of the continuum-subtracted peak [Call] 
flux. Therefore, even 20% changes in the continuum level 
result in only a ~3% change in the flux at any given wave¬ 
length. 


Additionally, we examined the consequences of phase 
on the velocity shift. Unfortunately, only one Ca-rich SN has 
sufficient data to examine any evolution. SN 2010et has two 
nebular spectra of reasonable quality that are separated by 
more than a week. At epochs of roughly 62 and 87 days after 
peak brightness, we measure a velocity shift of —330±60 and 
—430 ± 60 kms“^, respectively. These measurements differ 
by only 1.1 cr, but may represent a slight velocity gradient 
of —4 ± 3 kms“^day“^. 

SNe la have been s hown that their veloci ty shifts get 
redder with time (e.g., ISilverman et al.| 120131 '). eventually 
plateauing at a particular velocity as the ejecta become opti¬ 
cally thin. While this is a possibility for Ca-rich SNe, we note 
that SN 2012hn, which has the largest velocity offset and is 
blueshifted also ha s one of the latest spectra (-1-150 days; 
IValenti et al.ll20i3 '). 

Currently, there are not enough data to determine the 
velocity evolution of Ca-rich SNe with phase. Nonetheless, 
we note that the phase of a SN’s spectrum and its galactic 
offset are uncorrelated (correlation coefficient of r = 0.006), 
suggesting any potential correlation between offset and ve¬ 
locity is not caused by a correlation between phase and off¬ 
set. 


5.3 Measured Velocity Shifts 

In Figure [51 we present the continuum-subtracted [Call] 
AA7291, 7324 profiles for the full Ca-rich SN sample. There 
is a large diversity to the line profiles with roughly equal 
numbers having symmetric, skewed blue, and skewed red 
profiles. 

For each spectrum in Figure [51 we mark the emission- 
weighted velocity. We present these data in Table [AT] 

For f our SNe (2005E, 2010et, PTFOOdav, and 

PTFllbiil. IK12| measured [Call] velocity shifts. Their mea¬ 
surements deviate from ours by —80 kms“^ on average, 
with a median difference of —69 kms“^. The largest dif¬ 
ference is for PTF09da v, wh ere we measure a shift of 
—110 ± 60 kms“^, while I ki 3 measure 250 kms“^. While 
im do not describe their method for measuring the veloc¬ 
ity shift or list any uncertainties, this represents at most a 
5.9-a difference (using only the uncertainty associated with 
our measurement). All other differences are <200 kms“^ 
indicating a general agreement in these measurements. 

For 9 Ca-rich SNe, we could also measure an [O I] ve¬ 
locity shift, albeit with much larger uncertainties. Perform¬ 
ing a Bayesian Monte-Ca rlo linear re gression on the [Call] 
and [O I] velocity shifts dKellvl [2003) , we measure a slope 
and offset of 0.9 ± 0.9 and 0.2 ± 0.5, respectively. There¬ 
fore, velocity measurements from the [O I] feature are con¬ 
sistent with those of [Call], providing additional confidence 
that the [Call] velocity shifts measure any bulk offset to the 
ejecta. 

Five of 13 Ca-rich SNe have absolute velocity shifts of 
<300 km s“^ relative to their host galaxy rest frame, which 
are consistent with coming from galactic motion. However, 
8 SNe have larger shifts, up to 1700 km s“^ relative to their 
host galaxy rest frame. Such large velocity shifts must orig¬ 
inate from either an asymmetric SN explosion or extreme 
line-of-sight motion for the progenitor. 
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Figure 5. Left. Optical spectrum of SN 2003dg (black). Overplotted are template elliptical (red) and reddened {E{B — V) = 0.35 mag) Sc 
(blue) spectra. The galaxy spectra near [Call] AA7291, 7324 is not necessarily linear. Right. Corresponding continuum-subtracted [Call] 
AA7291, 7324 profile for SN 2003dg. The black, red, and blue spectra correspond to a linear continuum, elliptical galaxy, or reddened 
Sc galaxy subtracted from the SN spectrum. The resulting velocity shifts are plotted above the profile. For this example, the choice of 
continuum results in a systematic difference in the velocity shift of —240 and —80 kms“^for the elliptical and Sc templates, respectively. 


6 ANALYSIS 

6.1 The Velocity-Offset Correlation 

For the Ca-rich SN sample, there is a large range of ve¬ 
locity shifts (Section [5]). However, this range changes with 
projected offset. SNe with small physical projected offsets 
have a much larger range of velocity shifts than those with 
large projected offsets. In particular, within a projected off¬ 
set of 8 kpc, 4/7 Ca-rich SNe have absolute velocity shifts 
of >500 km s“^, while only 1/6 Ca-rich SN beyond 8 kpc 
has such a large absolute velocity shift (Figure 0. 

We also display these data as cumulative distribution 
functions (CDFs) for Ca-rich SNe within and outside a pro¬ 
jected offset of 8 kpc (Figure [SI which also shows the differ¬ 
ent line-of-sight velocity distributions for these subsamples. 
To examine the significance of this difference, we employ 
the Anderson-Darling statistical test. The Anderson-Darling 
test is similar to the Komogorov-Smirnov test except it is 
more sensitive to differences in the tails of distributions. 
Such a test is particularly useful for this analysis since we 
would naively expect any velocity shift distribution to be 
centered at roughly zero velocity and differences to be in the 
tails of the distributions. If the distributions are normal, the 
Anderson-Darling test reduces to the Komogorov-Smirnov 
test. 

The Anderson-Darling statistic for these two subsam¬ 
ples results in a p-value of 0.031, indicating that the Ca-rich 


SNe with small and large projected offsets are drawn from 
parent populations with different velocity shifts. This dis¬ 
tinction is relatively insensitive to the exact distance used 
to separate the subsamples. Any chosen separation ranging 
from 7.1 - 16.7 kpc produces significantly different subsam¬ 
ples. 

We also note that while we may be systematically miss¬ 
ing Ca-rich SNe near the centers of their host galaxies with 
redshifted velocity shifts (Section id.lll . the inclusion of these 
SNe, if they have a similarly extreme velocity distribution 
as those with blueshifted velocity shifts, would only make 
the tails of the distributions more discrepant and difference 
between the subsamples more significant. 

This result is incredibly indicative. First, if the veloc¬ 
ity shifts were primarily caused by the explosion or a bi¬ 
nary orbital velocity, there should be no correlation with 
projected offset. Therefore, the velocity shifts are almost 
certainly the result of extreme line-of-sight motion for the 
progenitor. To account for the SNe with extreme velocity 
offsets, some progenitor systems must be given a velocity 
kick of >1500 kms“^. 

There is an additional clue which indicates that the ve¬ 
locity shifts are not cau sed by the explosio n. Velocity shifts 
are seen in SNe la (e.g.. lMaeda et aDl2010l) and those shifts 
are thought to be a bulk offset of the core of the ejecta rela¬ 
tive to the outer layers. However, SN la explosions appear to 
be constrained to have relatively small velocity shifts com- 
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Projected Offset (kpc) Projected Offset (Isophotal Radii) 

Figure 7. Line-of-sight velocity shifts as a function of projected offset {left) and isophotal radial offset (right) for the Ca-rich sample. 
Velocity uncertainties are typically smaller than the data points. Only 4/13 SNe have a redshift, while only 1/6 SNe within two isophotal 
radii (and 8 kpc) have a redshift and 3/7 outside two isophotal radii (and 8 kpc) have redshifts. 


pared to the full-width at half-maximum (F WHM) of their 
emiss ion lines. Using a sample of 22 SNe la (jBlondin et al.l 
l2012h . we find that the most extreme velocity shifts are 
<20% that of their FWHM. For the Ca-rich SN sample, 
3/13 have more extre me velocity shif t s (rela tive to FWHM) 
than any SN la in the lBlondin et all (120121 ') sample. 

Figure |9] shows the CDFs of velocity shifts relative to 
FWHMs for the Ca-rich and SN la samples. The Anderson- 
Darling statistic for these two subsamples results in a p-value 
of 0.040, indicating that these samples are drawn from dif¬ 
ferent parent populations. This result implies that either 
Ca-rich SNe come from significantly more asymmetric ex¬ 
plosions than SNe la or there is an additional line-of-sight 
velocity component for some Ca-rich SNe. The latter is con¬ 
sistent with a correlation between velocity shift and pro¬ 
jected offset. 

The velocity-offset correlation further indicates that the 
line-of-sight motion for the progenitor must be correlated 
with the projected offset. The most obvious explanation is 
that most progenitor systems originate near the centers of 
their host galaxies, are given a large kick, and when they 
ultimately explode, have the velocity kick imprinted in the 
kinematics of the ejecta. For a progenitor that travels along 
our line of sight, there will be no projected offset for the SN, 
but its line-of-sight velocity offset will be large. Meanwhile, 
for a progenitor that travels perpendicular to our line of 
sight, the SN will have a large projected offset but a velocity 
offset close to zero (although the SN explosion or orbital 
motion of the progenitor system, as well as measurement 
errors may result in some velocity offset). 

Therefore, Ca-rich SNe progenitor systems do not typi¬ 
cally originate from globular clusters, in-falling dwarf galax¬ 
ies, the intracluster medium, or other large-offset compo¬ 
nents of a galaxy. 

From the velocity-offset data alone, it is possible that 


Ca-rich SN progenitor systems originate in the disk or halo, 
offset from the galaxy’s nucleus, and (in some cases) are 
kicked to large offsets. 

6.2 An Excess of Blueshifts 

The velocity offset distribution for the Ca-rich SN sample is 
skewed to the blue. Nine of the 13 SNe are blueshifted, while 
the remaining 4 are redshifted. This is particularly striking 
when examining the 6 SNe within two isophotal radii of 
their host galaxies (Figure (Tjl. Only a single SN in this sub¬ 
sample is redshifted. However, 3/7 of the Ca-rich SNe with 
larger offsets are redshifted. The likelihood of having such 
skewed distributions are only 17.5% and 18.8% for the full 
sample and small-offset subsample, respectively. This cal¬ 
culation accounts for the possibility that the distribution 
could be skewed in the opposite direction as well (if only 
considering blueshifts, the likelihood is half of that reported 
above). 

Although this is not a statistically significant result, it 
is noteworthy because of a physically motivated reason for 
such a difference. If the progenitor systems were kicked from 
the centers of their galaxies with velocities larger than the 
escape velocity of the galaxy, then the velocity shift will be 
blueshifted for SNe on the near side of the galaxy. In such a 
scenario, redshifted SNe will be on the far side of their host 
galaxies, and therefore will be more likely to suffer from dust 
extinction. Thus, Ca-rich SNe with redshifted velocity shifts 
will be harder to detect, resulting in under-representation in 
the full sample. 

This is further support for most Ca-rich SN progenitors 
originating from near the centers of their host galaxies. If 
their progenitors came from bound orbits (such as from a 
globular cluster population), there should be an equal num¬ 
ber of blueshifted and redshifted SNe on the near side of the 
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Figure 6. Continuum-subtracted [Call] AA7291, 7324 profiles of 
the Ca-rich SN sample on a velocity scale relative to the host- 
galaxy frame. The spectra are ordered by physical projected off¬ 
set. A circle indicates the emission-weighted velocity of each fea¬ 
ture. The uncertainty for each velocity measure is plotted, how¬ 
ever, most uncertainties are smaller than the plotted symbols. 
The offset and velocity (in 10^ km s“^) are listed next to each 
spectrum. 



Figure 8. Line-of-sight velocity shift CDFs for the Ca-rich SNe 
with a projected offset <8 kpc (black line) and >8 kpc (blue line). 
The SNe with small projected offsets have a wider distribution of 
line-of-sight velocity shifts (p = 0.031). 



Velocity Shift/FWHM 

Figure 9. CDFs of velocity shifts relative to nebular line FWHMs 
for the Ca-rich (black) and SN la (red) samples. Some Ca-rich SNe 
have much larger velocity shifts relative to the widths of their 
lines compared to the SN la sample. This indicates an additional 
velocity component. The Anderson-Darling statistic for these two 
subsamples results in a p-value of 0.040. 


galaxy, and therefore no noticeable difference in the observed 
population. Furthermore, for progenitors in-falling towards 
the host galaxy for the first time, we would expect more 
redshifted SNe, the opposite of what is observed. 
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Figure 10. Offset-angle CDF for Ca-rich SNe in edge-on galaxies 
(black). Also shown are theoretical distributions corresponding to 
a nuclear origin (blue line), a disk origin with the distance trav¬ 
eled by the progenitor system, d, relative to the major axis of the 
disk, r, of 0.11 or 0.4 (d/r; red and gold curves, respectively). The 
disk origin with dfr = 0.11 is inconsistent with the observed dis¬ 
tribution, while larger offsets relative to the height or the nuclear 
origin are consistent with the data. 


6.3 The Offset Angle Distribution 

A final observational indication of the origin of Ca-rich SNe 
from the SNe themselves is the angle between the SN and 
the disk of the galaxy, which we call the offset angle. A non¬ 
disk origin should have a flat offset angle distribution (all 
angles are equally likely). However, a disk origin will tend 
to have more SNe with small angles relative to the disk. 
The exact offset angle distribution for the disk depends on 
the average distance traveled by the progenitor system, d, 
relative to the major axis of the disk, r. For example, if 
d/r = 0, corresponding to the SNe exploding in the disk, 
then all offset angles would be zero. A larger value of d/r 
means that SNe will be further above the disk (on average) 
and thus the offset angle distribution will be flatter. 

In Figure [TOl we plot the offset angle distribution for the 
Ca-rich SNe with close to edge-on disks (axis ratios of <0.5). 
This corresponds to 5 SNe: SNe 2000ds, 2001co, 2003dg, 
2003dr, and 2005E. Of these 5, SNe 2000ds and 2003dr have 
offset angles >65°. The observed distribution is consistent 
with a flat distribution. 

To generate the theoretical disk distributions, we as¬ 
sumed an infinitely thin exponential disk perfectly edge on. 
For this stellar distribution, we chose random angles in three 
dimensions, and offset an amount d along those angles. We 
then measured the offset angle to the SN site. As expected, 
a smaller value of d/r results in a distribution more peaked 
to zero degrees. 

We use the Anderson-Darling statistic to determine 
what values of d/r are consistent with the data. Ca-rich 


SNe are inconsistent with coming from a disk population 
with ratios of d/r ^ 0.11. 

We cannot know the value of d for our SNe, however, 
we can measure the distance the SNe are above the disk, 
z, which must be a lower limit on d. For the edge-on sub¬ 
sample, there is a range of 0.04 ^ z/r 0.56. Therefore, a 
value of d/r « 0.4 may be appropriate, making a disk origin 
consistent with the observed angle offset distribution. 

More observations should further constrain any possible 
disk contribution. 


7 A BASIC MODEL 

Previous studies have shown that Ca-rich SNe have He in 
their ejecta (e.g.. lKawabata et al.ll2010l : IP10ll and come from 
older stellar populations (e.g., IPIOI) . These observations sug¬ 
gest an old progenitor system, and a massive-star progenitor 
is unlikely. The most obvious long-lived progenitor system 
with a large abundance of helium is one that contains a 
He WD or a hybrid He-C/O WD. As these stars are not ex¬ 
pected to explode in isolation, a favorable progenitor system 
is a He-C/O WD binary system. Simulations suggest that a 
merger (or surface detonation after mass transfer) for these 
syste ms can broadly rep r oduce the Ca-rich SN properties 


(e.g.. 

Kromer et al. 2010l: Shen et al. 201fll: IWaldman et al.l 

120111: 

Sim et al. 20121: Dessart & Hilliei] 2015|). 


Ca-rich SNe have an extended projected offset distri¬ 
bution (e.g., [Hi). While some have attempted to explain 
this as Ca-rich SNe having dwarf galaxy, globular cluster, 
and/or intracluster origins, these scenarios are inconsistent 
with other data. For instance, if there is a dwarf galaxy 
origin, there should not be such a strong correlation with 
nearby early-type galaxies. 

The strongest evidence against these previous scenarios 
is that there is no indication of any stellar populati on at 
the position of any Ca-rich SN to quite deep li mits 
IPerets et al.ll20fll : lK12l : iLvman et al.ll2013l . l2014|j . This last 
observation requires that the proge nitors of Ca-rich S Ne be 
born far from where they explode llLvman et al.llioTl '). 

Our kinematic data indicates that the progenitors of 
Ca-rich SNe originate near the centers of their galaxies. 
There is a strong correlation between the line-of-sight ve¬ 
locity shifts and projected offsets for the Ca-rich SN sample 
(Section [6T}. The measured velocity shifts for some SNe is 
much larger than typical stellar velocities and even larger 
than the escape velocity of a galaxy (Section [5]). In fact, 
the line-of-sight velocities of several Ca-rich SNe is larger 
than the expected m aximum velo city from a disk origin 
('^1000 km s“^; e.g.. iTauri^l2015h . but is consistent with 
that of Galactic hyper-vel ocity stars (HVSs; iBrown et al.l 
I 2 OI 4 I : [Palladino et al.|[20l3 l. From kinematic data alone, it 
is most likely that Ca-rich progenitor systems have been 
kicked after an interaction with a supermassive black hole 
(SMBH) at the center of their host galaxies. 

We have shown that almost all Ca-rich SN host galax¬ 
ies are merging, have recently merged, have characteristics 
indicative of a recent merger, and/or are in dense environ¬ 
ments where the probability of a merger is high (Section l3.1D . 
A binary SMBH can easily fill its loss cone, and therefore 
the HVS ejection rate of such systems are enhanced by a 
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factor of ^ 10 *^ iIYu fc Tremainel [200^ over those of single 
SMBHs re.g.. lHillslll9^ 

The ejection rate of hyper-velocity binary stars 
(HVBSs) is expected to be exceedingly low for single 
SMBHs, but perhaps half of all binary stars interacting with 
a binary S MBH are ejected at hyper velocitie s with the bi¬ 
nary intact et al.ll200^ : ISesana et al.ll2009lj . This mecha¬ 
nism is also expected to change the orbital properties of the 
binary. In particul ar, eccentricitie s of 0.9 - 1 can be excited 
in these systems (ILu et al ]|2^. Presumably, such orbits 
will again circularize, but in doing so, significantly reduce 
their separation. 

If a close encounter with a SMBH is necessary to pro¬ 
duce a Ca-rich SN progenitor system or to significantly re¬ 
duce its delay time, a recent galaxy merger should also sig¬ 
nificantly enhance the Ca-rich SN rate. In retrospect, one 
might have been able to predict that Ca-rich SNe should be 
preferentially hosted by galaxies that have undergone recent 
mergers. 

Additionally, we have shown that a large fraction of Ca- 
rich SN host galaxies also host an AGN (Section l3.2l) . After a 
merger, it is expected that many SMBHs will have their gas 
supplies replenished, and such galaxies will host an AGN. 
However, such activity can be relatively short-lived. The typ¬ 
ical duty cycle of an AGN is ~4 x 10^ years. Because AGN 
activity is relatively short-lived, any correlation between Ca- 
rich SNe and AGN activity suggests a physical connection. 
If one simply wants information to be transferred between 
the nucleus and SN location, that information must travel 
at a speed of >750 kms“^ for a separation of 30 kpc (sim¬ 
ilar to the separation between PTFllbij and the AGN in 
its host galaxy, 34 kpc). Notably, this velocity is larger than 
typical SN kicks. However, the expected hyper-velocity kicks 
from a SMBH interaction is typically, ~1500 kms“^, larger 
than, but within a factor of 2 of, the minimum speed for 
information travel in some Ca-rich galactic systems. 

With these strong pieces of evidence (helium in the sys¬ 
tem, an old stellar population, a lack of a stellar population 
at the position of the SNe, a correlation between line-of- 
sight velocity shifts and projected offset, a significant num¬ 
ber of host galaxies with a high-likelihood of recent mergers, 
and a high AGN fraction for the host galaxies), we suggest 
a possible progenitor system for Ca-rich SNe: a He-C/O 
WD binary system that has recently been ejected from its 
host galaxy after an encounter with a SMBH. The ejection 
rate will be significantly enhanced in galaxies with a binary 
SMBH, and therefore also in galaxies which have recently 
undergone a merger. Although it is possible for a Ca-rich 
SN to occur without an encounter with a SMBH, this must 
not be the dominant progenitor path for Ca-rich SNe. If Ca- 
rich SNe commonly came from systems that did not interact 
with a SMBH, we would not see the strong connections to 
mergers, AGN activity, or even the large offsets. Instead, 
the interaction with the SMBH must significantly enhance 
the SN rate, likely by significantly decreasing the delay time 
(perhaps from a typical delay time of greater than a Hubble 
time to one similar to the AGN duty cycle). 


7.1 Additional Model Constraints 

Such a progenitor model has other observational conse¬ 
quences consistent with current data. For instance, we would 



Host Galaxy Major Axis (kpc) 


Figure 11. Magnitude of velocity shifts as a function of host- 
galaxy size (as measured by the major axis) for the Ca-rich sam¬ 
ple. 


expect a flat offset angle distribution. Such a distribution 
is consistent with our observations (Section 16.311 . Addition¬ 
ally, as the SNe would be on unbound orbits, the SNe pro¬ 
jected on top of their galaxies would have the highest veloc¬ 
ity shifts. However, if there is dust reddening in the galaxy, 
we would expect redshifted SNe to be more difflcult to de¬ 
tect, resulting in an excess of SNe with blueshifted velocity 
shifts. Such an excess has also been observed (Section 16.21) . 

A final testable prediction of our model is that there 
should be an anti-correlation between the magnitude of the 
velocity shift and the size of the galaxy. Larger galaxies with 
larger gravitational potential wells should reduce the veloc¬ 
ity shift more than smaller galaxies. Alternatively, if the SN 
progenitors originated in satellite galaxies or globular clus¬ 
ters, the SN systemic velocity relative to their host galaxy 
should increase with increasing host-galaxy mass. 

Figure [TT] shows the magnitude of the velocity shifts as 
a function of host galaxy size for our sample. There is an ob¬ 
vious anti-correlation with the highest-velocity SNe coming 
from the most compact galaxies. While this is not a strong 
correlation (r = —0.29), the size of a galaxy is not a perfect 
proxy for gravitational potential, and cluster environments 
must further augment that potential, this is further obser¬ 
vational support for our model and evidence against a dwarf 
galaxy or globular cluster origin. 


7.2 Model Implications and Predictions 

Our model requires that the inherent Ca-rich SN rate be 
extremely low and an encounter with a SMBH must some¬ 
how substantially increase the rate. This is perhaps through 
changing the orbital dynamics of the system or by trigger¬ 
ing a physical change in a star in the system. There is a 
clear example of the former possibility. If the typical pro¬ 
genitor system is a binary system with a separation large 
enough such that a merger through gravitational radiation 
is longer than a Hubble time, then the inherent Ca-rich SN 
rate will be necessarily low. However, an interaction with 
a SMBH (perhaps preceded with other interactions in the 
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nuclear star cluster) could then reduce the merger time to 
significantly less than a Hubble time. 

We can estimate the post-kick merger time for these 
systems. If we assume that a typical kick velocity is 
1000 kms“^, a star system would have traveled for 30 or 
150 Myr to travel 30 or 150 kpc, respectively. While it is 
possible that some kicks result in bound orbits or that some 
SNe with large projected offsets also have relatively large 
line-of-sight offsets, these scenarios cannot be a significant 
portion of the population. Therefore, in our model, the typi¬ 
cal post-kick delay time is about 10-100 Myr. This timescale 
is similar to that of the AGN duty cycle. 

However, if this delay time is equivalent to a binary 
merger time, it must require a change in orbital separation 
by a factor of 10-100. This requirement may be reduced 
some if the binary is initially hardened by interactions in 
the nuclear star cluster. If this is the case, we would expect 
many more Ca-rich SNe in galactic nuclei. Although there 
are known systematics for detecting nuclear SNe, the lack 
of a substantial population of nuclear events in the current 
sample is a hurdle for such an explanation. 

Detailed calculations of the orbital parameters of poten¬ 
tial progenitor systems after ejection are required to measure 
a change in final separation and merger time. 

Another implication of the model is the relative rate of 
Ca-rich SNe in galaxies with and without recent mergers. 
Assuming that the progenitor systems are formed before a 
merger and all galaxies have a similar number of progenitor 
systems near their SMBHs, the relative Ca-rich SN rate per 
unit mass should be equivalent to the relative ejection rate, 
or ~10'*. Since the ratio of recent mergers to non-merging 
systems is roughly 10“^, the observed relative rate should 
be ~10^. That is, for every Ca-rich SN in a non-merging 
system, there should be 100 Ca-rich SNe in merging sys¬ 
tems. The current Ca-rich sample appears to have ~1/13 
non-merging systems is consistent with this estimate. How¬ 
ever, this estimate is only approximate and could be off by 
an order of magnitude. Similarly, the current sample may 
be highly biased or we may have incorrectly estimated the 
number of merging systems. This question should be revis¬ 
ited with more detailed calculations and a larger SN sample. 

Nonetheless, the relative rate of Ca-rich SNe to SNe la 
has been measured to be ~5%. If we assume that Ca-rich 
SNe occur exclusively in merging systems, the Ca-rich SN 
rate must be ~5 times that of the SN la rate in these galax¬ 
ies. The SN la rate in SO galaxies is roughly 0.1 SNe la 
per century llLi et al. llmll) or 10“® year"^ A 6.5 x 10® Mq 
SMBH with a full loss cone has an ejection rate of ~2 stars 
per year. Therefore, Ca-rich SN progenitor systems must 
represent ~0.25% of all ejected systems in these galaxies. 
Considering the potentially high surv ival rate of co mpact 
binary systems with binary SMBHs (ILu et al. lli^ , this 
order-of-magnitude calculation appears to have the Ca-rich 
SN rate compatible with the ejection rate. 


8 DISCUSSION AND CONCLUSIONS 

Ca-rich SNe are rare and peculiar low-luminosity Type I SNe 
with He in their spectra near maximum brightness and rela¬ 
tively strong Ca emission at late times. While many of these 
SNe are within one isophotal radius of their host galaxies, a 


significant percentage occur far from their host galaxies — 
some farther than any other detected SN. 

The Ca-rich SNe with small projected offsets tend to 
have much larger line-of-sight velocity shifts as determine 
from the emission-weighted velocity of forbidden lines. The 
correlation between projected offset and line-of-sight veloc¬ 
ity shifts indicates that the progenitors of Ca-rich SNe have 
been kicked from near the centers of their galaxies. 

The highest velocity shifts are too large for a disk ori¬ 
gin, but are consistent with a kick from an interaction with 
a SMBH. A nuclear origin is further supported by a slight 
excess of blueshifted velocity shifts with small offsets, an off¬ 
set angle distribution consistent with isotropic kicks from a 
galactic nucleus, and an anti-correlation between the mag¬ 
nitude of the line-of-sight velocity shifts and galaxy size. 

The host galaxy population has a significant number of 
merging/disturbed galaxies, SO galaxies, and brightest clus¬ 
ter/group galaxies — all of which are more likely to have had 
recent mergers than typical galaxies. Such an observation 
must be linked to the prevalence of Ca-rich SN progenitor 
systems. A recent galaxy merger enhancing the interaction 
rate between Ca-rich SN progenitor systems and a SMBH 
is one possible reason for this observation. Such a scenario 
is expected to boost the SMBH ejection rate by a factor of 
>^10^, making observed and theoretical rates broadly consis¬ 
tent. 

The fraction of Ca-rich SN host galaxies with AGN is 
high. The AGN fraction is yet another indication of recent 
mergers. The timescale for AGN activity also places a soft 
upper limit on the delay time since galaxy merger of -^4 x 
10 ^ years and an independent lower limit on the velocity of 
ejected systems. 

The host-galaxy morphology dist ribution is indicative 
of an older progenitor population (IPlOh . The He in the early- 
time spectra of Ca-rich SNe suggests a significant amount 
of He is in the progenitor system. As the progenitor system 
is likely old, the source is likely a He WD. 

As a possible progenitor hypothesis, we propose the fol¬ 
lowing: the progenitor system for a Ca-rich SN is a binary 
system consisting of a He WD with a C/O WD. This sys¬ 
tem must have a typical merger time longer than a Hubble 
time, which is most likely achieved with a low-mass dou¬ 
ble WD system. This system is kicked to high velocity after 
interaction with a SMBH in the center of its host galaxy. 
During this interaction, the binary is significantly hardened 
such that the delay time is reduced to -^50 Myr. Such a 
possibility seems plausible for binary systems ejected after 
interacting with a binary SMBH, whic h is also expect ed to 
significantly reduce the merger time dLu et ahllioOTll . Af¬ 
ter traveling, on average, several to tens of kpc, the system 
mer ges a nd results in a Ca-rich SN. 

IPld examined the possibility of the progenitor of 
SN 2005E being ejected as a HVS, however they only exam¬ 
ined this for a potential high-mass (M > 8Mq) progenitor. 
They found that such a progenitor system was unlikely, but 
did not address the possibility of a low-mass binary progen¬ 
itor ejected at hyper velocity. 

If this scenario is correct, the progenitor systems of Ca- 
rich SNe represent the first indication of extragalactic HVSs 
and the first detection of HVBSs ejected from a galactic 
center. These peculiar SNe would provide the best way to 
observe HVSs beyond the local group. Ca-rich SNe may be 
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an excellent tracer of the SMBH population in general and 
binary SMBHs in particular. Ca-rich SN host galaxies may 
provide an excellent input catalog for future binary SMBH 
studies. Moreover, the Ca-rich SN rate may be an indirect 
tracer of the galaxy merger rate. 

We caution that the details of our toy model are cur¬ 
rently poorly constrained. The current sample size is small 
and likely biased against nuclear events (and perhaps against 
SNe with the largest projected offsets). An order of magni¬ 
tude larger sample, which can be achieved in the next few 
years with upcoming surveys combined with dedicated spec¬ 
troscopic resources, will provide necessary kinematic con¬ 
straints to infer the underlying properties of the Ca-rich SN 
progenitor population. 
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APPENDIX A: OBSERVATIONS OE PTFllkmb 

PTEllkmb was discovered PTF on 24.24 August 2011 and 
was spectroscopically classifie d as a SN Ib from a spectrum 
obtained on 28 August 2011 (iGal-Yam et aLlboilf) . 

PTFllkmb was observed on 28.52 August 2011 (PI Fil¬ 
ippenko; Program U048LA) and 26.37 November 2011 (PI 
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Figure A12. Optical spectra of PTFllkmb (black curves) from 
28 August 2011 (top) and 26.37 November 2011 (bottom). Com¬ 
parison spectra of SN 2005 E fro m 15.35 January 2005 (top) and 
11.26 March 2005 Ibottom: IPlOll are displayed in red. All spectra 
have been corrected for Milky Way reddening. The continuum is 
slightly different at early times, but this is likely the result of being 
at somewhat different phases or having different temperatures. 


Kulkarni; Program G219L A) with the Low Resolution Imag¬ 
ing Spectrometer ILRIS: loke et aDll995^ mounted on the 
Keck I telescope. The August spectrum had a 450 s integra¬ 
tion, while the November spectrum had 1200 s blue-channel 
and 2x570 s red-channel integrations. The August spectrum 
was originally used to classify PTFllkmb as a SN Ib. 

We obtained these data through the Keck Observatory 
Archive. Standard CCD processing and spectrum extraction 
were accomplished with IRA^ The data were extracted us¬ 
ing the optimal algorithm ofE orn 3 (ll986l) . Low-order poly¬ 
nomial fits to calibration-lamp spectra were used to establish 
the wavelength scale, and small adjustments derived from 
night-sky lines in the object frames were applied. We em¬ 
ployed our own IDL routines to flux calibrate the data and 
remove telluric lines using the well-exposed cont i nua of spec- 
trophotometric standards (|Wade &: Horii^l 19881 : iFolev et al.l 
I2OO3I : ISilverman et al.|l2012l) . The spectra are presented in 
Figure IA12I 

As the comparisons to SN 2005E in Figure IA12I show, 
PTFllkmb is clearly a member of the Ca-rich class of SNe. 
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Table Al. Ca-rich SN Host-galaxy Demographics 


SN 

Host 

Morphology 

AGN 

z 

Scale 

(kpc/") 

Major 
Axis {") 

Axis Ratio 
(fe/a) 

Galactic 
RA. (°) 

2000ds 

NGC 2768 

E 

LINER 

0.0045 

0.113 

65.70 

0.460 

-87.5 

2001CO 

NGC 5559 

Sb (BGG) 

BPT Comp 

0.0172 

0.383 

18.55 

0.300 

60 

2003H 

NGC 2207 

Sbc (interacting) 

AGN 

0.0091 

0.171 

36.25 

0.680 

70 

2003dg 

UGC 6934 

Scd 


0.0183 

0.405 

11.15 

0.320 

-35 

2003dr 

NGC 5714 

Scd (BGG) 


0.0075 

0.184 

28.30 

0.160 

80 

2005E 

NGC 1032 

so/a 

Radio Excess 

0.0090 

0.173 

33.80 

0.400 

70 

2005CZ 

NGC 4589 

E (BGG) 

LINER 

0.0066 

0.158 

34.75 

0.750 

92.5 

2007ke 

NGC 1129 

E (BCG; merger) 


0.0173 

0.340 

38.15 

0.880 

70 

2010et 

CGCG 170-011 

E (group) 

BPT AGN 

0.0233 

0.495 

7.00 

0.800 

-30 

2012hn 

NGC 2272 

SO (BGG) 


0.0071 

0.142 

25.50 

0.740 

-55 

PTF09dav 

2MASS J22465295-I-2138221 

Disturbed Spiral 


0.0371 

0.735 

4.35 

0.740 

-85 

PTFllbij 

IC 3956 

so? (BGG) 

AGN 

0.0347 

0.720 

5.45 

0.900 

-20 

PTFllkmb 

NGC 7265 

SO (BGG) 


0.0170 

0.343 

26.25 

0.760 

-10 


Table A2. Ca-rich SN Sample Demographics 


SN 

SN Offset 
(") 

SN Offset 
(kpc) 

SN Offset 
(Isophotal Radii) 

Angle 
Offset (°) 

Velocity Shift 
(km s“^) 

2000ds 

33.4 

3.77 

1.08 

76.67 

-190 (30) 

2001CO 

18.5 

7.07 

1.22 

13.03 

-470 (530) 

2003H 

51.0 

8.73 

1.52 

22.47 

-180 (130) 

2003dg 

4.1 

1.66 

0.39 

6.05 

-1120 (40) 

2003dr 

14.4 

2.65 

2.89 

65.09 

560 (40) 

2005E 

140.3 

24.27 

6.07 

27.78 

-140 (30) 

2005CZ 

13.4 

2.12 

0.41 

24.07 

-1030 (140) 

2007ke 

49.2 

16.71 

1.31 

17.31 

470 (120) 

2010et 

76.0 

37.64 

10.95 

9.44 

-430 (70) 

2012hn 

47.4 

6.73 

2.48 

75.50 

-1730 (70) 

PTF09dav 

56.8 

41.75 

16.09 

52.29 

-110 (60) 

PTFllbij 

47.8 

34.42 

8.80 

10.00 

580 (100) 

PTFllkmb 

437.5 

150.05 

20.46 

56.36 

260 (20) 









